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Abstract 

We revisit the solvability of quadratic semimartingale BSDEs and propose a Lipschitz- 
quadratic regularization procedure. In the first step, we prove an existence and uniqueness 
result for Lipschitz-quadratic BSDEs. The proof is inspired by but excluded from [16]. 
A corresponding stability result and a Lipschitz-quadratic regularization are developed to 
solve quadratic BSDEs. The advantage of our approach is that much weaker conditions 
guarantee the existence and uniqueness results. 
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1 Preliminaries 


In this paper, we are concerned with the solvability of K.-valued backward stochastic differ¬ 
ential equations (BSDEs) driven by continuous local martingales which take the form 


Y t=Z + f t (1 T d(M) s F(s,Y s ,Z s )+g s d{N) s ) - £ ( Z s dM s +dN s ), 


(1) 


where M and N are strongly orthogonal continuous local martingales. We are particularly 
interested in the above equations with quadratic growth, i.e., the generator F is quadratic 
in Z and g is not identical to 0. 

BSDEs of this type have been intensively applied to mathematical finance and stochastic 
control; see Mania and Schweizer [12], or Hu et al [6] in Brownian setting. In its theoretical 
aspect, Karoui and Huang [9] obtains the solvability with Lipschitz-continuous generators. 

Later, Tevzadze [16] studies the existence and uniqueness of a bounded solution, by assum¬ 
ing quadratic growth and local Lipschitz-continuity. Morlais [14] extends the stability-type 
argument in Kobylanski [10] to quadratic BSDEs driven by continuous local martingales. 

Based on this work, Mocha and Westray [13] proves existence and uniqueness results with 
convex generators and exponential moments integrability. 

A close inspection of this line of study, however, reveals that their assumptions are 
quire demanding. For example, the stability-type argument in Morlais [14] can be used 
only if the BSDE is not quadratic in N, i.e., g. = 0. When g is a constant process, an 
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exponential transform can be used to kill the quadratic term g ■ (N ). But one has to 
sacrifice the flexibility of the generators, especially for unbounded solutions. For this point, 
the interested readers shall refer to [9], [14], [13]. When g is a bounded process, some results 
are obtained by Tevzdaze [16], but rather restrictive. For example, existence results are 
obtained only for particular quadratic generators, and equations with Lipschitz-continuous 
generators are not studied. 

Having understood these literature and their drawbacks, we develop a Lipschitz-quadratic 
regularization technique to answer the question of existence and uniqueness under more flex¬ 
ible assumptions. In the first step, we study BSDEs with Lipschitz-continuous generators 
and quadratic growth in N, by adapting the fixed point arguments in Tevzadze [16]. These 
equations, due to this particular structure, are called Lipschitz-quadratic. Viewing this 
result as a basic building block, we then derive a corresponding monotone stability result 
to faciliate our study of more general quadratic BSDEs. The regularization therein is called 
Lipschitz-quadratic, as contrary to the Lipschitz regularization in [14], [13]. It turns out 
that all the results, including existence, uniqueness and stability results of bounded and 
unbounded solutions can be obtained with weaker conditions. 

This paper is organized as follows. In Section 2, we prove an existence and uniqueness 
result for Lipschitz-quadratic BSDEs. Based on this result, we establish a monotone stabil¬ 
ity theorem in Section 3. As a byproduct, the existence of a bounded solution is immediate. 
In Section 4, we study existence, uniqueness and stability results for unbounded solutions, 
using a localization procedure. Finally, Section 5 reviews the change of measure result 
studied in Mocha and Westray [13]. 

Let us close this section by introducing all required notations. We fix the time horizon 
0 < T < Too, and work on a filtered probability space (f2, J 7 , (J 7 t)te[o,71>I l> ) satisfying the 
usual conditions of right-continuity and P-completeness. _Fo is the P-completion of the trivial 
cr-algebra. Any measurability will refer to the filtration (J~t)te[o.T] ■ I 11 particular, Prog 
denotes the progressive cr-algebra on 12 x [0, T] . We assume the filtration is continuous , in the 
sense that all local martingales have P-a.s. continuous sample paths. M = (M 1 ,..., M d ) T 
stands for a fixed d-dimensional continuous local martingale. By continuous semimartingale 
setting we mean: M doesn’t have to be a Brownian motion; the filtration is not necessarily 
generated by M which is usually seen as the main source of randomness. Hence in various 
concrete situations there may be a continuous local martingale strongly orthogonal to M, 
which we denote, as in (1), by N. 

Here we clarify all notions in (1). We set 1 := (1,..., 1) T . £ is an K-valued JV-measurable 
random variable, F : D x [0,T] x K x I d -r is a Prog®B(R) ® £>(R d )-measurable 
random function and g is an R-valued Prog-measurable bounded process. J 0 (Z s dM s +dN s ), 
sometimes denoted by Z ■ M + N, refers to the vector stochastic integral; see Shiryaev and 
Cherny [15]. The equations defined in this way encode the matrix-valued process (M) 
which is not amenable to analysis. Therefore we rewrite the BSDEs by factorizing (M). 
This procedure separates the matrix property from its nature as a measure. It can also be 
regarded as a reduction of dimensionality. 

There are many ways to factorize (M); see, e.g., Section III. 4a, Jacod and Shiryaev 
[8]. We can and choose A := arctan (JV =1 (Af*)). By Kunita-Watanabe inequality, we 
deduce the absolute continuity of with respect to A. Note that such choice 

makes A continuous, increasing and bounded. Moreover, by Radon-Nikodym theorem and 
Cholesky decomposition, there exists a matrix-valued Prog-measurable process A such that 
(M) = (A T A) ■ A. As will be seen later, our results don’t rely on the specific choice of A but 
only on its boundedness. In particular, if M is a d-dimensional Brownian motion, we may 
choose A t = t and A to be the identity matrix. 
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The second advantage of factorizing (M) is that 


1 T d(M) s F(s,Y s ,Z s ) 


l T AjA s F{s,Y s ,Z s )dA s , 


where f(t,y,z) := l T AjA s F(s, y, z) is R-valued. Such reduction of dimensionality makes 
it easier to formulate the difference of two equations as frequently appears in comparison 
theorem and uniqueness. Hence, we may reformulate the BSDEs as follows. 

BSDEs: Definition and Solutions. Let A be an R-valued continuous nondecreas¬ 
ing bounded adapted process such that (M) = (A T A) • A for some matrix-valued Prog- 
measurable process A, / : Q x [0, T] x R x ! d —» R a Prog<g>B(R) <8> B(R d )-measurable 
random function, g an R-valued Prog-measurable bounded process and £ an R-valued Tt~ 
measurable random variable. The semimartingale BSDEs are written as 


F t =£+ 


J (/(s, Y s , Z s )dA s + g s d(N) s s ) — j ( Z s dM s + dN s ). 


( 2 ) 


We call a process (Y, Z,N) or (Y, Z-AI + N) a solution of (2) , if Y is an R-valued continuous 
adapted process, Z is an Revalued Prog-measurable process and N is an R-valued contin¬ 
uous local martingale strongly orthogonal to M, such that P-a.s. J Q T Zjd(M) s Z s < +oo 
and / q T \f(s,Y s ,Z a )\dA s < +oo, and (2) holds P-a.s. for all t £ [0,T], 

Note that the factorization of (M) gives f Q Zjd(M) s Z s = f Q \\ s Z s \ 2 dA s P-a.s. Hence 
we don’t distinguish these two integrals in all situations, Zjd(M) s Z s < +oo P-a.s. 
ensures that Z is integrable with respect to M in the sense of vector stochastic integration. 
As a result, Z ■ M is a continuous local martingale. M and N being continuous and strongly 
orthogonal implies that (M l ,N). = 0 for i = 1 ,...,d. We call / the generator , £ the terminal 
value and (£, f Q T | f(s, 0, 0) \dA s ) the data. In our study, the integrability property of the data 
determines the estimates for a solution. The conditions imposed on the generator are called 
the structure conditions. For notational convenience, we sometimes write (/, g,£) instead 
of (2) to denote the above BSDE. Finally, (2) is called quadratic if / has at most quadratic 
growth in z or g is not indistinguishable from 0. 

To finalize, we introduce the rest notations which will be used throughout this paper. 
<C stands for the strong order of nondecreasing processes, stating that the difference is 
nondecreasing. For any random variable or process Y, we say Y has some property if this is 
true except on a P-null subset of H. Hence we omit “P-a.s” in situations without ambiguity. 
Define sgn(x) = For any random variable X, define HATH^ to be its essential 

supremum. For any cadlag adapted process Y, set Y a j. '■= Yt — Y s and Y* := sup tg r 0 T ] |Yi|. 
For any Prog-measurable process H , set \H\ s ,t '■= /* H u dA u and \H\ t := T stands 

for the set of all stopping times valued in [0,T] and S denotes the space of continuous 
adapted processes. For later use we specify the following spaces under P. 

• S°°: the space of bounded processes Y £ S with ||y|| := ||F’*|| 00 ; S°° is a Banach 
space; 

• A4: the set of continuous local martingales starting from 0; for any Revalued Prog- 
measurable process Z with J Q T Zjd(M) s Z s < +oo, Z ■ M £ M\ 

• Ai p (p > 1): the set of M £ M with 

\\M\\ mp := (E[(M)|])' <+oo; 

in particular, M 2 is a Hilbert space; 
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• M bmo : the set of BMO martingales M £ M with 

\\M\\ bmo := sup||E[(M) r , T |JVp|| ; 
reT 

M.bm° is a Banach space. 

M 2 being a Hilbert space is crucial to proving convergence of the martingale parts in 
the monotone stability result of quadratic BSDEs; see, e.g., Kobylanski [10], Briand and 
Hu [2], Morlais [14] or Section 3. Other spaces are also Banach under suitable norms; we 
will not present these facts in more detail since they are not involved in our study. 

Finally, for any local martingale M, we call {cr, l } ngN + cTa localizing sequence if a n 
increases stationarily to T as n goes to +oo and M. A(Tn is a martingale for any n £ N + . 

2 Bounded Solutions of Lipschitz-quadratic BSDEs 

This section takes one step in solving quadratic BSDEs and consists in the study of equations 
with Lipschitz-continuous generators. In contrast to El Karoui and Huang [5], we allow the 
presence of g • (TV). We point out that similar results for linear-quadratic generators have 
been studied by Tevzadze [16], but the case of Lipschitz-continuity is not available in that 
work. Due to its importance for regularizations of quadratic BSDEs, we study existence 
and uniqueness results for equations of this particular type in the first step. To this end, 
we assume 

Assumption (A.l) There exist /3,y > 0 such that + || 11/(-, 0,0)| | T || < +oo and 

/ is Lipschitz-continuous in {y, z), i.e., P-a.s. for any t £ [0,7), y,y' £ R, z, z' £ R d , 

I f(t,y,z) - f{t,y',z')\ < 0\y - y'\ + 7 |A t (z - z')\. 

Due to the presence of g • (TV), we call the BSDE (f,g,£) satisfying (A.l) Lipschitz- 
quadratic. Given (A.l), we are about to construct a solution in the space 98 := S°° xA4 BMO 
equipped with the norm 

\\(Y, Z ■ M + N)\\ := (\\Y\\ 2 + \\Z.M + N\\ 2 bmo )\ 

for (y, Z ■ M + TV) £ S°° x M bmo . Clearly {38, ||-||) is Banach. As a preliminary result, 
we claim that the existence result holds given sufficiently small data. 

Theorem 1 (Existence (i)) If (f,g,£) satisfies (A.l) with 

liaL + 8|||l/(-.0.0)|| T || 2 oo <iexp(-||A||(8/3 2 ||A||+8 7 2 )) (3) 

and P-a.s. \g.\ < g := then there exists a solution in {38, ||-||). 

Proof. To overcome the difficulty arising from the Lipschitz-continuity, we use Banach 
fixed point theorem underran equivalent norm. Set p > 0 to be determined later. For any 
X £ L°°,y £ S°° and M £ M bmo , set ||-X'|| 00># , := ||e** r .X’|| 00 , ||y|| p := \\eS A Y\\ and 

ll M ll BMO.p '■= \\e* A - M Wbmo'i for {Y,Z ■ M + N) £38, set 

IKE Z • M + TV)|| p := (||y|| 2 + \\Z ■ M + N\\ 2 BMO p ) *. 
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Since A is bounded, ||-|| is equivalent to the original norm for each space. Hence {38, ||-|| ) 
is also a Banach space. For any R > 0, define 


:= {{Y,Z ■ M + N) £ 38 : \\{Y,Z ■ M + N)\\ p < R}. 

We show by Banach fixed point theorem that there exists a unique solution in B/j with 
R = To this end, we define F : (B#, ||-|| ) —>• {38, ||-|| ) such that for any {y, z ■ M +n) £ 
Bfi, (Y, Z ■ M + N) := F {{y, z ■ M + n)) solves 


Yt= £ + J t (f( s ’y^ z s)dA s + g s d(n)s) - ( Z s dM s + dN s ). 


Indeed, such (Y, Z, N) uniquely exists due to martingale representation theorem. Moreover, 
by standard estimates, {Y, Z ■ M + N) £ {38, ||-|| ). 

(i). We show F(B#) C Hr. For any t £T, Ito’s formula applied to e pA Y 2 yields 


e pAr \Y T \ 2 + pE 
<ll?llL,p + 2E 


e pAa YgdA s 


Tr 


+ E 




(zJd(M) s Z s + d{ N), 


T t 


e pA “ |Y s ||/(s, y s , z s )\dA s 


T t 


•2E 


e^|Y s ||. 9s |d(n), 


T t 


(4) 


By (A.l), 


\ Y a\\f{s,y a ,z a )\ < |F s ||/(s,0,0)1 +P\Y s \\y s \ + 7 |Y 5 ||A s z s |. 

We plug this inequality into (4) and estimate each term on the right-hand side. Using 
2 ab < \a 2 + 8 b 2 gives 


2E 


e pAa \Y s \\f{s,0,0)\dA s 


T r 


<|nip+8E 


<^mip+8|||i/(-.o,o)i 


e^ As l/( s , 0,0)|eL4 s 
2 


t t 


1 2 


T II oo,p’ 


2/?e[ / e pA °\Y s \\y s \dA s T r ] < -\\y\f + 8/3 2 e[ / e? A ‘\Y s \dA 


T t 


<^ll 2 /||p + 8/3 2 P||E[ / e pAs \Y s \ 2 dA s 


Tr 


2 7 e[ / e pA *\Y s \\\ s z s \dA s T t ] <-\\z • M\\% MOtf) + 8 7 2 e[ / e pA ‘\Y a \HA. 


T r 


2E 


e^|Y s || 5s |(7V), 


T t 


<^m„+8 9 2 E 


e^ A ‘d{N) t 


T t 



\\Y\\ 2 p + 8~g 2 \\n\\ 


4 

BMO,p m 


5 









































Set p := 8/3 2 ||H|| + 87 s so as to eliminate E[ J^ e pAs Y^dA s \ J>] on both sides. Hence (4) 
gives 


oP A t 


\Yt 


-E 


oP A ‘ 


( Zjd(M) s Z s +d(N ) ; 




< IICH 2 oo.p + 8|||l/(-,o,o)||. 


Vii? 


>,p 1 '-'i i| x IIoo,p ' 4 II" !I P 

II 2 ' M\\bmo, p ) + 8 g 2 ||n|| BMOjp . 


lip 1 ll~ ‘‘ , ‘ t \\tfMU,pJ 1 II * v WUMU,p' (5) 

Taking essential supremum and supremum over all t£ T, and using the inequality 


1 


z-m + n)\\; < \\y\\; V ||Z• M + N\\ BMOf . 


< sup 
rer 


?P A T 


ini 


-E 


0 pA 


(z s T d(M ) 6 




T-r 


we deduce by transferring ^||Y|| p to the left-hand side of (5) that 

||(Y, Z ■ M + N)\\l < 4||e||L,p + 32||||/(.,0,0)|| r ||^ p + ±(\\y\\ 2 p + \\z ■ M\\ 2 bmo J + 32g 2 ||n 

< 4iienL, P + 32|| 11/(-, 0, ojiyiL., + \r 2 +32 g 2 R 4 . 

Thanks to (3), g = | and R = we verify from the above estimate that 


\\(Y,Z,N)\\ p <R. 

(ii). We prove F : (B#, ||-|| p ) —> (B#, ||-|| p ) is a contraction mapping. By (i), for i = 1,2 
and any ( y\z l ■ M + ri 1 ) £ B^, we have ( Y l , Z l ■ M + N l ) := F ((y l , z l ■ M + ri)) £ B fi . 
For notational convenience we set Sy := y 1 — y 2 and 6z,Sn,S(n),5Y,5Z,SN,S(N), etc. 
analogously. By the deductions in (i) with minor modifications, we obtain 


l -\\(5Y,5Z-M + 5N)\\ 2 p < 1 - 


\\^-M\\ 2 bmo J + j\\6Y\\ 


+ 4g 2 sup 
reT 


E 


RA 
e 2 - 


d\S(n) s 


T-r 


1 2 


( 6 ) 


Kunita-Watanabe inequality and Cauchy-Schwartz inequality used to the last term gives 


E 


£ A 

O -'T-f 


d\5(n) s 




< E 


3 d(Sn) 


Tr 

in 2 


E 


e : ‘ Aa d(n 1 + n 2 ) s 


Rr 


^ ll^ n ll BMO,p ' 2(ll^ 1 II BMO,p + ll n2 |l BMO,p) 


— Il^ n ll BMO,p ‘ ^R ? 

where the last inequality is due to || (y z , z 1 • M + ri) || < R, i = 1,2. Hence ( 6 ) gives 


\\(6Y,6Z.M + 6N)\\l<±'"'" 2 


< 


||i5z • M\\ BMO p ) + 645 2 i? 2 ||(5n|| BMO p 

(i + 64ff 2 i? 2 ) || (Sy, Sz ■ M + Sn )|| p 


< ^|| (Sy,Sz-M + Sn )|| p , 


4 

BMO,p 
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i.e., F : (B fl , ||-|| p ) -s> ( Bflj IML) is a contraction mapping. The existence of a solution 
in Bjj thus follows immediately from Banach fixed point theorem. Finally, since ||-|| is 
equivalent to ||-|| for 38, the solution also belongs to (38, ||-||). 


□ 


From now on we denote (38, ||-||) by 38 when there is no ambiguity. In the spirit of 
Tevzadze [16], we extend this existence result so as to allow any bounded data. To this end, 
for any Q equivalent to P we define iS°°(Q) analogously to S°° but under Q. This notation 
also applies to other spaces. 

Theorem 2 (Existence (ii)) If (f,g,f) satisfies (A.l), then there exists a solution of 

(f,g,€) in 38. 

Proof, (i). We first show that it is equivalent to prove the existence result given \g.\ < ^ 
P-a.s. Suppose that g is bounded by a positive constant g , that is, |g.| < g P-a.s. Observe 
that, for any 9 > 0, (Y, Z, N ) is a solution of (/, g, £) if and only if (9Y, 9Z, 9N) is a solution 
of (f e ,g/9,9tf), where f 9 (t,y,z) := 9f(t, |, |). Obviously f e verifies (A.l) with the same 
Lipschitz coefficients as /. If we set 6 := 8 g, then \g./6\ < | P-a.s. and hence satisfies the 
parametrization in Theorem 1 (existence (i)). Therefore, we can and do assume \g. | < | 
P-a.s. without loss of generality. 

(ii). Since + || 11/(•, 0,0)|| T || < +oo, we can find n G N + such that 

n n 

/(*. 0 , 0 )=$:/*(*, 0 , 0 ), 

i= 1 i =1 


where, for each i < n, is a J^-measurable random variable, /* : O x [0, T] x ffi. x K. d —> K. 
is ProgOS(K.) <g> S(R d )-measurable and 

rilL+8|||ir(-,0,0)|| T £<^exp(-||A||(8/3 2 ||A||+8 7 2 )). 

Set f'(t, y , z) := f(t, y, z) — f(t, 0, 0) and (F°, Z° ■ M + N°) G 38 such that ||(F°, Z° ■ M + 
iV°)|| = 0. Now we use a recursion argument in the following way for i = 1,..., n. 

By Theorem 1, there exists a solution (Y l , Z l ■ M + N l ) G 38( Q l ) to the BSDE 


2 — 1 2-1 


Yi = c+ (ns, o, o)+ f( S , J2 y s j ,J2 z s ) - /'(«> E Y i, E z d) 


dA, 


j—0 j =o 


J=0 j= 0 


ffs d(lV l ) s - / (ZldM s + dNl), 


where 


dP 


:=£ 


2-1 

(2 g-J2 NJ 

j =0 


Note that the equivalent change of measure holds due to the fact that N 3 G M. bmo for 
j < i — 1 and Theorem 2.3, Kazamaki [9]. By Girsanov transformation and Theorem 3.6, 
Kazamaki [9], N l := N l + 2g ■ (N l , AT- 7 ) and Z r ■ M belong to M. BMO . This further 
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implies (N*) = (W) and W = N i + 2g-(N i , N j ). Hence (Y\ Z^M + N*) £ S8 solves 

/ T i i i—1 i —1 

(r ( s , o, o)+ ns, y yLY z i ) - /'(«> E E z ii) dA ° 

3 =0 j=0 3=0 3=0 

+ [ g s d((N') s + 2(N i .,Y NJ )s) - [ ( ZldM s +dNl ). 

Jt v j=Q ' Jt 

Hence a recursion argument gives (Y 1 , Z\ N l ) for i = 1, n. 

Define Y := £”=0 Y\ Z ■= EEo and N : = EEo Clearly (Y, Z • M + iV) G 9S. 
We show (' Y,Z,N) solves (/,g,£). I n view of the definition of /', we sum up the above 
BSDEs to obtain 

Y t =Z + Jt ((f(s,0,0) + f'(s,Y a ,Z a ))dAs+g a d{N) a ')- J (6Z s dM s + dSN s ). 

To conlcude the proof we use f'(s, Y s , Z s ) := f(s, Y s , Z s ) — f(s, 0,0). 


□ 


We continue to show that comparison theorem and hence uniqueness also hold given 
Lipschitz-continuity. Similar results in different settings can be found, e.g., in [12], [ 6 ], [14], 
[16]. 

Theorem 3 (Comparison) Let (Y, Z ■ M + N), (Y', Z' ■ M + N') £ S°° x M bmo be 
solutions of {f,g,£), (/ , ,g ,/ ,£ / )> respectively. IfP-a.s. for any ( t,y,z) £ [0,T] x R x R d , 
f{t,y,z ) < f'(t, y, z), g t < g' t , (<£' and (f,g,f) verifies (A.l), then F-cl.s. Y. < Yf. 

Proof. Set 5Y := Y — Y' and 5Z,5N,6(N),S £, etc. analogously. For any r £ T, P-a.s. 
/ < /' and g. < g[ imply by ltd’s formula that 


SYtAr =SYr+ [ (/(a , Y„ Z s ) - f'(a, Yf, Z’ s ))dA s + [ g s d(N) s - [ g' s d(N') s 

JtAr JtAr JtAr 

- I ( SZ s dM s + dSN s ) 

J tAr 

< SY T + [ (/(a, Y„ Z s ) - /(a, Yf, Z' s ))dA s + [ g' s dS{N) s - [ ( SZ s dM s + dSN s ) 

J tAr J tAr J tAr 

= 8Y t + f (/3 S 6Y S + (Xs'js) 1 ^ {^ s SZ s ))dA s + f g' s d6(N) s — f (5Z s dM s + dSN s ), 

JtAr JtAr JtAr 


(7) 


where /3 (M-valued) and 7 (Revalued) are defined by 

f(s,Y s ,Z s )-f(s,Yf,Z s ) 


Ps ■= I{5F S ^0} ; 

7s '■ = \\ s 5Z a ^ 0} 


SY s 


(f(s,Yf,Z a )-f(s,Yf,Z' B ))6Z s 

\XsSZ s \ 2 


and 0 := (0,...,0) T . Note that 7 can be seen as defined in terms of discrete gradient. By 
(A.l), j3. and f Q jJd(M) s j s are bounded processes, hence 7 • M £ A4 bmo . Given these 




facts we use a change of measure to attain the comparison result. To this end, we define a 
BMO martingale 

A-.= ' Y -M + g' -{N + N'). 

In view of Theorem 2.3 and Theorem 3.6, Kamazaki [9], we define 



Hence SN - g’ ■ S(N) and SZ ■ M — ( 7 T A T A SZ) ■ A belong to M bmo (Q). 
and P-a.s. < 0 give 


SY t < JG 


/3 s 5Y s dA s 


Tt 


< 


p s SY s dA a 




Therefore, (7) 


Hence we obtain by Gronwall’s lemma that P-a.s. SY t < 0. Finally by the continuity of Y 
and Y', we conclude that P-a.s. Y. < Yf. 


□ 

As a byproduct, we obtain the following existence and uniqueness result. 

Corollary 4 (Uniqueness) //(/, g,f) satisfies (A.l), then there exists a unique solution 
in 3). 

Proof. This is immediate from Theorem 2 (existence (ii)) and Theorem 3 (comparison 
theorem). 


□ 

3 Monotone Stability and Bounded Solutions of Quadratic BSDEs 

In this section, we prove a general monotone stability result for quadratic BSDEs. Let 
us recall that Morlais [14] uses a stability-type argument for the existence result after 
performing an exponential transform which eliminates g ■ (N) . But a direct general stability 
result is not studied. Our work fills this gap. 

Secondly, as a byproduct of the stability property, we construct a bounded solution via 
regularization through Lipschitz-quadratic BSDEs studied in Section 3. This procedure 
is also called Lipschitz-quadratic regularization in the following context. Note that our 
definition of “Lipschitz-quadratic” is different from those in [16], [1]. To begin our proof, 
we give the assumptions for the whole section. 

Assumption (A.2) There exist /? > 0, 7 > 0, an K + -valued Prog-measurable process a 
and a continuous nondecreasing function ip : R + —» M+ with y>( 0 ) = 0 such that + 

IIHalloo < +°° an d P-a.s. 

(i) for any t G [0, T], (y, z) 1 — > f{t , y, z ) is continuous; 

(ii) / is monotonic at y = 0, i.e., for any (t, y , z) £ [0, T] x R x R d , 

sgn {y)f{t,y,z) < a t +a t /3\y\ + ||A t z| 2 ; 
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(iii) for any (t, y, z) G [0, T] x 1 x l d , 

\f(t,V,z)\ < a t + a t <p(\y\) + ^\X t z\ 2 . 

We continue as before to call (£, |a|r) the data. (A.2)(ii) allows one to get rid of the 
linear growth in y which is required by Kobylanski [10] and Morlais [14]. Assumption of 
this type for quadratic framework is motivated by Briand and Hu [2]. Secondly, our results 
don’t rely on the specific choice of tp. Hence the growth condition in y can be arbitrary as 
long as (A.2) (i) (ii) hold. 

Given (A.2), we first prove an a priori estimate. In order to treat (Z ■ M) and g ■ (N) 
more easily, we assume P-a.s. \g. \ < 7 for the rest of this paper. 

Lemma 5 (A Priori Estimate) If (/, g , £) satisfies (A.2) and (Y, Z ■ M + N) G S°° x M 
is a solution of (f,g,£), then 

PI < ||e /3WT (|£| + 


and 


||Z ■ M + A || BMO < Cb, 


where Cb is a constant only depending on ||£|| , IHcxItIP • 

Proof. Set u(x ) := exp ^ 73 ^)~ 1 ~' 1 ' J . The following auxiliary results will be useful: u(x) > 
0,u'(x) > 0 and u"(x ) > 1 for x > 0; u(| • |) G C 2 (K) and u"{x) = 7 u’(x) + 1. For any 
r, er G T, Ito’s formula yields 

u(\Y TAa \) =u(\Y a \) + [ u'(|F s |)sgn(F s )dF s -i [ u"{\Y s \)(zJd(M) s Z s + d(N) s ). 

Jt/\ct " JrAcr ^ ' 

By (A.2) (ii), 

sg n(Y s )f(s,Y s ,Z s ) < a s + a s /3\Y s \ + ^\\ S Z S \ 2 . 

Note that iu'(|F s |) - |u"(|y s |) = -±, ff s u'(|F s |) - \u"{\Y s \) < and «'(|y a |) < 

Hence, using these facts to the above equality yields 


1 

2 



( 'zJd(M) s Z s + d(N ) s ) 


0 7||X|| 


< — 

7 


u , (|y s |)(a s + a s P\Y s \)dA s 


u'(\Y s \)sgn(Y s )(Z s dM s +dN s ). 


To eliminate the local martingale, we replace a by its localizing sequence and use Fatou’s 
lemma to the left-hand side. Since Y* and |o|t are bounded random variables, the right- 
hand side has a uniform constant upper bound. Hence, we have 

1 e 7pl e7 Pl 

-K[{Z.M + N) T , T \T T \<— + —{l+mi)\\W\ T \\ 00 - ( 8 ) 


Now we turn to the estimate for Y. We fix s G [0, T\ and for t G [s, T], set 


H t 


:= exp 
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We claim that H is a submartingale. By Tanaka’s formula, 


d\Y t \ = sgn (Y t )(Z t dM t + dN t ) - sgn (Y t )(f(t, Y t ,Z t )dA t + g t d(N ) t ) + dL°(Y), 
where L°(Y) is the local time of Y at 0. Hence, Ito’s formula yields 

dH t = jHte 0 ^ [sgn{Y t )(Z t dM t + dN t ) 

+ ( - sgn {Y t )f(t, Y t , Z t ) + at + atP\Y t \ + | e fl“l-*|A t Z t | 2 )dA t 
+ ( - sgn (Y t )g t + l e 0 ^)d(N) t + dL° t (Y) . 

By (A.2)(ii) and \g.\<^ again, [Ht)te\s,T] is a bounded submartingale. Hence, 


|F S | < -lnE[fr T |^]. 


Thanks to the boundedness, we have 

m < pw-dei + HrjiL- 

Finally we come back to (8) and obtain the estimate for Z ■ M + N. 


□ 

Given the norm bound in Lemma 5, we turn to the main result of this section: the 
monotone stability result. Later, as an immediate application, we prove an existence result 
for quadratic BSDEs by Lipschitz-quadratic regularization. To start, we recall that .Ad 2 
equipped with the norm ||iV4"||_ /V12 := E[(M)t] 2 for M G At 2 is a Hilbert space. 

Theorem 6 (Monotone Stability) Let(f n ,g n ,C) n &^ satisfy (A. 2) associated with (a, /?, 7 , <p), 
and ( Y n , Z n ■ M + N n ) be their solutions in 3d, respectively. Assume 

(i) Y n is monotonic in n and ff 1 — £ — > 0 P-a.s. with sup n ||^ ra || 00 < +00; 

(ii) P-a.s. for any t G [0 ,T], gf - g t —> 0; 

(iii) P-a-s- for any t G [0, T] and y n —> y, z n —> f n (t, y n , z n ) —> f(t, y, z ). 

Then there exists a process iY, Z ■ M + N ) G SS such that Y n converges to Y P-a.s. 
uniformly on [0, T) and ( Z n ■ M + N n ) converges to (Z ■ M + N) in M 2 as n goes to + 00 . 
Moreover, ( Y,Z,N ) solves (f,g,£). 

Proof. Without loss of generality we only consider Y n to be increasing in n. By Lemma 
5 (a priori estimate), 


sup||Fl + sup||Z n • M + N n \\ BMO < c b , (9) 

n n 

where q, is a constant only depending on /?, 7 ,sup n ||£"ll 00 ) IIMrlloo- We rely intensively on 
the boundedness result in (9) to derive the limit. 

(i). We prove the convergence of the solution sequences. Due to (9), there exists a 
bounded monotone limit Yt := lim„ Yf l , a subsequence indexed by {nfe}fc £ N+ Q N + and 
Z ■ M + N G A4 2 such that Z nk ■ M + N nk converges weakly in M 2 to Z ■ M + N as k goes 
to Too. The remaining task is to show Z ■ M + N is the At 2 -limit of the whole sequence. To 
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this end, we define u(x) := exp ^ 87 ^ 7 a 871 1 ■ Recall that u(x) > 0 ,u'{x) > 0 and u"{x) > 0 
for x > 0; u G C 2 (R) and u"{x) = 87 u'(x) + 1. For any m G {nfc}fceN+ ,n G N + with m > n, 
define 5Y m ’ n := Y m - Y n , 5Y n :=Y -Y n and 5Z m ’ n , 6Z n , SN m ’ n , <5TV”, etc. analogously. 
By Ito’s formula, we have 


E [u(5Y 0 m ’ n )\ -E[u(<5C’")] =E / u\5Yr n ){r{s,Yr,ZT)-f n {sX:\Z:))dA s 


L J 0 


+ E 


JO 




-gy(N n )s) 

r T 


L Jo 


u"(SY s m ’ n ) (5Z ™’ n ) 1 d(M) s (5Z™ ,n ) + diSN™' 71 ), 


( 10 ) 


Since / m and / n verify (A.2) associated with (a, ft, 7 , </?), we have 

!/"*(«, 17 *, Z 7 *)-/"(«, y/*,z a ")| 

<a'+||A s Z 7| 2 + ||A s ^| 2 

< + ^{\X s 6Z™’ n \ 2 + \X S 5Z™\ 2 + \X S Z S \ 2 ) +1 (\X S 6Z:\ 2 + \X S Z S \ 2 ) 

<a' s + ‘?l\X s 6Z?’ n \ 2 + ^(|A S ^"| 2 + \X S Z S \ 2 ), 


where 

a' s ■= 2a s (l + ip(c b )) > 2a s + a s g>{\Y™\) + a s < / ?(|Y s m |). 

Moreover, 

g m d(N m ) - g n d(N n ) < |d(TV m ) + | d(N n ) 

< ^ d(6N m ’ n ) + ^ (d(SN n ) + d(N)). 


Plugging the above inequalities into (10), we deduce that 


E 


L Jo 


1 „ 37 , 

-u -— W 

2 2 


y8Y s m ' n )\X s 6Z™’ n \ 2 dA s 
57 


+ E 


L Jo 


1 // 37 , 

-u ——w 

2 2 


) (6Y™’ n )d(5N m ’ n ) z 


< E[u(^ m ’ n )] +E[ / u'{8Y™’ n )(a' s + -^(|A S <^ S ”| 2 + |A s ^ s | 2 ))dA. 


+ E 


L «/0 


« , (^ Tn,n )Y(d<^)a+d(JV)a) 


( 11 ) 


Due to the weak convergence result and convexity of z 1 —> \z\ 2 , N 1 —> (TV), we obtain 


E 


£ (l u "-^ u, )(SY s n )\x t z :\ 2dA s 

' £ ik u " -Y u ') {SY ° n)d{6Nn)s 


L Jo 
E 


< liminfE 


< liminfE 


£ Q u " _ ^ u 'Y 5Y r n )\^ tZ 7' n \ 2d A. 

J o ( 8Y s m ’ n )d(8N m ’ n ) 
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We then come back to (11) and send m to +oo along {nfc}fcgN+. Taking the above inequal¬ 
ities into account and using u'(8Y™' n ) < u'(8Y s n ) to the right-hand side, (11) becomes 


E 


X Q u " - Y u ') {6Y ^ XsZ °\ 2dAs 


-E 


l u " 


[ SA 

<«[»«')] +E [j u'(SY, n )(a', + ^.(\X,SZH 2 + \\,Z,\ 2 )) dA - 


52e 

2 


u'(SY s n )(d(8N n ) s + d(N) s ) 


( 12 ) 


Since u"{x) — 87 u'(x) = 1, rearranging terms give 


Ie[(6N2) 2 ] + ±e[ [ \\ s 5Z:\ 2 dA s 


'-J 0 


<E [u(5C)]+E[J^ u\5Y s n )(^a , s + ^\X s Zs\ 2 )dA s \+^E[^ u\SY s n )d(N) s ]. (13) 


Finally, by sending n to +00 and dominated convergence we deduce the convergence. 

(ii). We prove (Y. Z ■ M + N) G 3% and solves Here we rely on the same 

arguments as in Kobylanski [10] or Morlais [14] and omit the details here. In addition to 
their deductions, we need to prove the u.c.p convergence of g n ■ (N n ), which holds if 


lim E 

n—)• 00 


/ ( gy(N n )s 

Jo 


g s d(N) s ) 


*1 


= 0 . 


Indeed, by Kunita-Watanabe inequality and Cauchy-Schwartz inequality, 


E 


{g n s d{N n )s - gsd(N) s ) 


= E 


(gy((N n ) s - (N) s ) + ((£ - g s )d(N) s ) 


< ^E[(N n - N) T ] 2 E[(N n + N) t ] 2 +E 


(g?-g s )d(N) t 


< 7CftE[(AI n — N)t] 5 + E / Iff" — gs\d(N) 


We then conclude by A4 2 -convergence of N n and dominated convergence used to the second 
term. Finally Z ■ M + N G J\4 BMO by Lemma 5 (a priori estimate). 

For decreasing Y n , we take m G N + ,n G {nfc}fe g pj+ with n > m and conclude with 
exactly the same arguments. 


□ 

There are several major improvements compared to existing monotone stability results. 
First of all, in contrast to Kobylanski [10] and Morlais [14], we get rid of linear growth in y 
by merely assuming (A.2), and allow g to be any bounded process. Secondly, we treat the 
convergence in a more direct and general way than Morlais [14]. 

Another advantage concerns the existence result. Thanks to Section 2 and Theorem 
6, we are able to perform directly a Lipschitz-quadratic regularization without exponential 
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transforms; this is in contrast to Morlais [14]. One can also benefit from our stability result 
in obtaining the existence results for unbounded solutions with more flexible assumptions; 
see Section 4. 

Proposition 7 (Existence) If (f,g,f) satisfy (A. 2), then there exists a solution in 98. 

Proof. We use a double approximation procedure and use Theorem 6 (monotone stability) 
to take the limit. Define 

f n ' k (t,y,z) : = inf {f + (t, y', z') + n\y - y'\ + n\\ t {z - z')l} 

y',z' 

- inf {f~(t,y',z') + k\y-y'\+k\\t(z-z')\}. 

y',z' 

By Lepeltier and San Martin [11], f n,k is Lipschitz-continuous in (■ y,z ); as k goes to +oo, 
fn,k conver g es increasingly uniformly on compact sets to a limit denoted by / n, °°; as n 
goes to Too, f n ’°° converges increasingly uniformly on compact sets to /. 

By Corollary 4, there exists a unique solution ( Y n,k , Z n ’ k -M + N n ’ k ) £ 98 to ( f n ’ k ,g , £); 
by Theorem 3 (comparison theorem), Y n,k is increasing in n and decreasing in k, and is 
uniformly bounded due to Lemma 5 (a priori estimate). We then fix n and use Theorem 6 
to the sequence indexed by k to obtain a solution ( Y n ,Z n ■ M T N n ) £98 to (/"’°°, <?, £)• 
Due to the P-a.s. uniform convergence of Y n,k we can pass the comparison property to Y n . 
We use Theorem 6 again to conclude. 


□ 

Remark. In contrast to Kobylanski [10], the existence of a maximal or minimal solution is 
not available (yet) given (A.l) as the double approximation procedure makes the comparison 
between solutions impossible. 

There is also a rich literature on the uniqueness of a bounded solution of quadratic 
BSDEs; see, e.g., [10], [12], [6], [14]. Roughly speaking, they essentially rely a type of 
locally Lipschitz-continuity and use a change of measure analogously to Section 2. The 
proof in our setting is exactly the same and hence omitted to save pages. 

To end this section, we briefly present various structure conditions used in different 
situations. 

Assumption (A.2') There exist /3 > 0,7 > 0, an K + -valued Prog-measurable process a, 
and a continuous nondecreasing function <p : R + —> R + with 99 ( 0 ) = 0 such that P-a.s. 

(i) for any t £ [0, T], ( y , z) 1 — > f{t, y, z) is continuous; 

(ii) / is monotonic at y = 0, i.e., for any (f, y, z) e[0,T]xIx R d , 

sgn (y)f(t,y,z) < a t +P\y \T ^\X t z\ 2 ; 

(iii) for any (t, y, z ) £ [0, T] x R x R d , 

\f{t,y,z)\ <a t + ip(\y\) T ^|A t 2 | 2 . 
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Given bounded data, (A. 2') implies (A.2). Indeed, 

sgn {y)f(t,y,z) < a t V 1 + (a t V 1)/%| + ||A t 2| 2 , 

\f(t,y,z)\ < a t V 1 + (a t V l)<p{\y\) + ||A t z| 2 - 

Hence (A.2') verifies (A.2) associated with (a V1, /3, 7 , </?). However, given unbounded data, 
(A.2') appears to be more natural and convenient. This will be discussed in detail in Section 
4. 

In particular situations where the estimate for f Q \f(s, Y s , Z s )\dA s is needed, e.g., in 
analysis of measure change (see Section 5), there has to be a linear growth in y, which 
corresponds to the following assumption 

Assumption (A. 2") There exist f > 0, 7 > 0, an Revalued Prog-measurable process a 
such that P-a.s. 

(i) for any t G [0, T], (y, z ) 1 — > f(t , y, z ) is continuous; 

(ii) for any (t, y, z ) G [0, T] x R x R d , 

\f(t,y,z)\ < a t +P\y\ + ^|A t z| 2 . 

Indeed, (A.2") enables one to obtain the estimate for / Q T |/(s, Y s , Z s )\dA s via 
\f(s,Y s ,Z s )\dA s < \a\ T +m\\Y* + 1(Z-M) t . 

4 Unbounded Solutions of Quadratic BSDEs 

This section extends Section 2, 3 to unbounded solutions. We prove an existence result 
and later show that the uniqueness holds given convexity assumption as an additional 
requirement. We point out that similar results have been obtained by Mocha and Westray 
[13], but our results rely on much fewer assumptions and are more natural. Analogously 
to section 3, we give an a priori estimate in the first step. We keep in mind that P-a.s. 

| g. | < ^ throughout our study. 

Lemma 8 (A priori estimate) If satisfies (A. 2') and ( Y,Z ■ M + N) G S x M 

is a solution of (/, g , £) such that the process 



is of class T>, then 



|U|<ilnE exp ^ 7 e /3As ’' r |£| + 7 J e 0As ' u a u dA 1 




Proof. We fix s G [0,T], and for t G [s,T], set 


H t := exp (^)e^ As ' t \Y t \ + 7 f e^ Ae ' u a u dA u ^ . 


(14) 


(15) 
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We claim that H is a local submartingale. Indeed, by Tanaka’s formula 


d\Y t \ =sgn (Y t )(Z t dM t + dN t ) - sgn (Y t )(f(t,Y t , Z t )dA t + g t d{N) t ) +dL° t (Y), 
where L°(Y) is the local time of Y at 0. Hence, Ito’s formula yields 


dH t = 'yH t e > 


fi As 


sgn (Y t )(Z t dM t + dN t ) 


+ ( - sgn (Y t )f(t,Y t , Z t ) + a t + fi\Y t \ + ^e 0A ^\X t Z t \ 2 )dA t 
+ ( - sgn(y t )<?t + \e^)d(N) t + dL° t (Y) . 

By (A.2')(ii), H is a local submartingale. To eliminate the local martingale part, we replace 
t by its localizing sequence on [s,T], denoted by {T„} ngN +. Therefore, 

\Y S \< -lnE[U TArn |j r s ] 


< — InE 

7 


fTAr n 


exp ^ 7 e /3As ’ TAr ” |Ftat„| +7 J 
Finally by class T> property we conclude by sending n to +oo. 


e P A *,u a dA, 


Ts 


□ 


We then know from Lemma 8 that exponential moments integrability on |£| + \cx\t is a 
natural requirement for the existence result. 


Remark. (A.2') addresses the issue of integrability better than (A.2). To show this, let us 
assume (A.2). We then deduce from Lemma 5 and corresponding class V property that 


\Y S \ < - InE 

7 


exp + 7 / e ^H a ’ u a u dA l 


A 


(16) 


Obviously, in (16), even exponential moments integrability is not sufficient to ensure the 
well-posedness of the a priori estimate. For more dicusssions on the choice of structure 
conditions, the reader shall refer to Mocha and Westray [13]. 


Motivated by the above discussions, we prove an existence result given (A. 2') and ex¬ 
ponential moments integrability. Analogously to Theorem 7, we use a Lipschitz-quadratic 
regularization and take the limit by the monotone stability result in Section 3. The a priori 
bound for Y obtained in Lemma 8 is also crucial to the construction of an unbounded 
solution. 

Theorem 9 (Existence) If (f, 9,0 satisfies (A. 2 ') and e /3AT (|^| + | oi\t) has exponential 
moment of order 7 , i.e., 


E 


exp (7e /3A ' r (|£| + Mt)) 


< + 00 , 


then there exists a solution verifying (14). 
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Proof. We introduce the notations used throughout the proof. Define the process 


X t := — InE 

7 


exp^e^dCI + Mr)) 


Ft 


Obviously X is continuous by the continuity of the filtration. For m,n £ N + , set 

T m ■= inf {t > 0 : \a\ t + X t >m} AT, 
a n := inf {t > 0 : \a\t >n}AT. 

It then follows from the continuity of X and |a|. that r m and a n increase stationarily to T 
as m,n goes to +oo, respectively. To apply a double approximation procedure, we define 

f n ' k (t, y, z ) := I {t<a n } inf {f + (t, y , z') + n\y-y'\+ n\X t (z - z') |} 

- y',z' 

- I{t< CTfc } inf {f~(t,y', z') + k\y - y'\ + k\X t (z — ^ 7 )|}, 


and £" ,fc := A n — A k. 

Before proceeding to the proof we give some useful facts. By Lepeltier and San Martin 
[11], f n ’ k is Lipschitz-continuous in (y,z)\ as k goes to +oo, f n ’ k converges decreasingly 
uniformly on compact sets to a limit denoted by /"’°°; as n goes to +oo, / n, °° converges in¬ 
creasingly uniformly on compact sets to F. Moreover, | |/ ra,fc (-, 0,0)| | T and £ n,k are bounded. 

Hence, by Corollary 4, there exists a unique solution (Y n ’ k , Z n,k ■ M + N n,k ) € 38 to 
(f n ' k ,g,£, n ' k )\ by Theorem 3 (comparison theorem), Y n,k is increasing in n and decreasing 
in k. Analogously to Proposition 7, we wish to take the limit by Theorem 6 (monotone 
stability). 

However, \f n ’ k {-, 0, 0)|t and £ n,k are not uniformly bounded in general. To overcome 
this difficulty, we use Lemma 8 (a priori estimate) and work on random interval where Y n,k 
and \f n ’ k (-, 0,0)|. are uniformly bounded. This is the motivation to introduce X and r m . 
To be more precise, the localization procedure is as follows. 

Note that (/"’ fc , g, £ n ’ fc ) verifies (A.2') associated with (a, f3, 7 , ip). Y n,k being bounded 
implies that it is of class V. Hence from Lemma 8 we have 


7 

< — In E 

7 

<X t . 


\Y t n ’ k | < ilnE[exp( 7 e^|C fe H-7y e 0A '-a s l {s < an ^ k} dA a 


Ft 


exp (^/e /3At ’ T \t;\ + j J e^ At ' T a s dA^j 


Ft 


(17) 


In view of the definition of r m , we have 

\Y t n *J<X tATm <m, 

\\r k (-,0,0)\\ Tm <\I [0 , rm] a\ Tm <rn. (18) 

Hence | |/"’ fe (-, 0, 0)11 and Y n,k are uniformly bounded on [0, r m ]. Secondly, given ( Y n,k , Z n,k - 
M + N n ’ k ) which solves (/"’ fc , g^ n ’ k ), it is immediate that (Y%? m , (Z n ’ k ■ M + N n ’ k ). ATm ) 
solves (I[ 0 ^ Tm }(t)f n ’ k (t,y,z),g,YJY k ). We then use Theorem 6 as in Proposition 7 to con¬ 
struct a pair ( Y m , ( Z m ■ M + N m )) which solves (/, g , sup„ inffc YjY k ), i.e., 

^ f‘Tm ^ _ _ f‘Fm 

Y t m = sup inf Y™^ k + / (F(s,Y s m ,Z™)dA s + g s (N m ) s ) - / (Z?dM, + dN„). 

( 19 ) 
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Moreover, Y m is the P-a.s. uniform limit of Y"r m and Z m ■ M + N m is the _A/( 2 -limit of 
( Z n,k ■ M + N n ’ k )./\ Tm as k,n go to + 00 . Hence 


Y%£ = Y™ Tm P-a.s., 
I{ t < Tm }A t Z™ +1 = A t Z™ dA ® dP-a.e, 

P-a.s. 

Define (Y, Z, N ) on [0, T] by 

Yt := I {t<r m }Yt + E WurJ YT, 

m> 2 

Zt ■= I {t<T m }Zt + E 1 

m> 2 

N t -.= ItK^+EWi,^’ 

m>2 


( 20 ) 


By (20), we have H Arm = X ATm ) an d ^V-Ar m = ^.Ar m - Hence we 

can rewrite (19) as 


Y tATrn 


pTm 

sup inf Y™^ + / 

n k m Jtf\T m 


(f(s, Y s , Z s )dA s + g s d(N ) s ) 



(. Z s dM s + dN s ). 


By sending to to + 00 , we prove that (Y, Y, A/) solves (/, g, ^). By (17), we have 


|It | = | sup inf Y) n,fe | < — InE exp [^e^ At ' T |^| + 7 / e^ At ' s a s dA t 
n k 7 L V 


•Ft 


□ 

Compared to Mocha and Westray [13], we prove the existence result under rather milder 
structure conditions. For example, (A.2')(ii) gets rid of linear growth in y and allows g to 
be any bounded process, which has been seen repeatedly throughout this paper. Secondly, 
in contrast to their work, the assumption that dA t CAdt , where ca is a positive constant, 
is not needed. Finally, they use a regularization procedure through quadratic BSDEs with 
bounded data. Hence, more demanding structure conditions are imposed to ensure that 
the comparison theorem holds. On the contrary, the Lipschitz-quadratic regularization is 
more direct and essentially merely relies on (A.2') which is the most general assumption to 
our knowledge. For the differences, the interested reader shall refer to [14], [13]. 

Due to the same reason as in Proposition 7, the existence of a maximal or minimal 
solution is not available. 

Remark. Analogously to Hu and Schweizer [7], one may easily extend the existence result 
to infinite-horizon case. In abstract terms, given exponential moments integrability on 
exp(/3A 00 )|o!| 00 , we regularize through Lipschitz-quadratic BSDEs with increasing horizons 
and null terminal value. Using a localization procedure and the monotone stability result 
as in Theorem 9, we obtain a solution which solves the infinite-horizon BSDE. 

As a result from Lemma 8 , we derive the estimates for the local martingale part. To 
save pages we only consider the following extremal case. 
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Corollary 10 (Estimate) Let (A. 2') hold for (/, g , £) and e^ Ar (|£| + |q;|t) has exponen¬ 
tial moments of all orders. Then any solution (Y, Z , N) verifying (14) satisfies: Y has 
exponential moments of all order and Z • M + N G Ai p for all p > 1. More precisely, for 
all p > 1 , 



exp (pje 0AT (\£| + Mr)) 


and for all p > 1, 


E 


pT p 

( J ' (zJd{M) s Z s + d(JV> s )) T ] < cE[exp (4p 7 e^(|£| + |a| T )) 


where c is a constant only depending on p, 7 . 


Proof. The proof is exactly the same as Corollary 4.2, Mocha and Westray [13] and hence 
omitted. 


□ 

Let us turn to the uniqueness result. We modify Mocha and Westray [13] to allow g to 
be any bounded process rather than merely a constant. A convexity assumption is imposed 
so as to use 0-technique which proves to be convenient to treat quadratic terms. We start 
from comparison theorem and then move to uniqueness and stability result. Similar results 
can be found in Briand and Hu [2] for Brownian setting or Da Lio and Ley [3] from the point 
of view of PDEs. To this end, the following structure conditions on (f,g,£) are needed. 

Assumption (A.3) There exist /3 > 0,7 > 0 and an R + -valued Prog-measurable process 
a such that P-a.s. 

(i) for any t G [0, T], (y, z ) 1 —> f(t , y, z ) is continuous; 

(ii) / is Lipschitz-continuous in y, i.e., for any (t,z) G [0, T] x R d , y,y' G R, 

I f(t,y,z) - f(t,y',z )| < (3\y - y '|; 

(iii) for any (t, y) G [0, T\ x R, z 1 —> f(t , y, z) is convex; 

(iv) for any (t, y, z ) G [0, T] x R x R d , 

\f(t,y,z)\ < a t +/3\y\ + ||A t ^| 2 . 

We start our proof of comparison theorem by observing that (A.3) implies (A.2'). Hence 
existence is ensured given suitable integrability. Likewise, we keep in mind that P-a.s. 

Ifl’-I — 2 ‘ 

Theorem 11 (Comparison Theorem) Let (Y, Z ■ M + N), (Y ', Z' ■ M + N') G S x A4 be 

solutions of (/, g,^), (f , g ', ff), respectively, andY*,{Y')*, \cx\t have exponential moments 
of all orders. If P-a.s. for any (t,y,z) G [0,T] x R x R d , f (t,y,z) < f'(t,y,z), g t < g' t , 
g' t >0, £ <£' and {f,g,£) verifies (A.3), then P-a.s. Y. < Yf. 

Proof. We introduce the notations used throughout the proof. For any 6 G (0,1), define 

Sf t :=f(t,Y',Z' t )-f'(tX,Z^ 

SgY :=Y — 0Y’, 

SY :=Y- Y', 
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and SgZ, SZ , 6gN,6N, etc. analogously. Moreover, define 


Pt ■= I{<5„Y t #0} 


f(t,Y t ,Z t )-f(t, 9 Yf,Z t ) 

SeY t 


By (A.3)(ii), p is bounded by fd for any 0 £ (0,1). Hence \p\r < /3||A||. By Ito’s formula, 
e' p '‘SeYt = eM T 6 e Y T + J e^F e s dA s + J e W* ( g s d{N) s - 9g' s d(N') s ) 
eM'(6gZ a dM B + d6 e N a ), 

Jt 

where 

F e s = f(s, Y S ,Z S ) - 0f\s , y;, Z') - p a S g Y a , 

= 66fs + ( f(s , y 8 , z 8 ) - /(s, y;, z 8 )) + {f(s, y;, z 8 ) - a/( s , y;, zj)) - p 8 fl e y 8 . ( 21 ) 

We then use (A.3)(ii)(iii) to deduce that 

f(s, Y a , Z a ) - f(t, Y', Z a ) = f(s, Y S ,Z S ) - f(s, fly;, Z a ) + /(«, fly/, Z 8 ) - f(s, YJ, Z a ) 

= p s SgY s + f(t , fly;, z a ) - /(«, y;, z 8 ) 

< Ps 6gY s + (1 - fl)^|y;i, 

/(*, n', - of (a, y;, zj) = /( S , y;, az' + (1 - a)^|) - a/(*, y;, zj) 


<(i-A)/( s ,y;,^) 

< (i - 0K + (i - fl)^|y;i + 


7 


-|A s ^y s | 2 . 


2(1-A) 1 

We also note that IP-a.s. 6f s < 0. Hence plugging these inequalities into (21) gives 
F e s < (1 - fl)(a 8 + 2(3\Y'\) + ^-^—\\ s 6 e Z s \ 2 . 

We then perforin an exponential transform to eliminate both quadratic terms. Set 


( 22 ) 


c := 


7e- 


MM 


1-0 ’ 

Pt ■= exp (ce^SgYt). 


By Ito’s formula, 

Pt=P T +Jt cPj»\° (f 6 s - c ^l\S e Z s \ 2 yA s 


C 1 | , / np\p\a \ 

+ j cP a eW‘ [g s d(N) s - 9g' s d(N') s - —d(5 e N) s ) 
- f cP s e\ p \ s ( 5 gZ s dM s + d 5 gN s ). 
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For notatioiial convenience, we define 


G t ~ 


■■= cP t eM*(F t e -^\Z°\ 2 ), 


H t := / cP 


By (22), we have 


where 


J cP s e M ° (g s d(^V) s - 9g' s d(N') s - -d(N e ) s ). 

G t = cP t eW‘((l-0)(a t + 2/3|F/|)) < P t J t , 


J* := 7 e 2/3 'l A ll(a t + 2/3|F t '|). 

We claim that H can also be eliminated. Indeed, 

d(S e N) = d(N) + d 2 d(N') - 29d(N, N') 

» d(N) + 9 2 d{N') - 9d{N) - 9d(N') 
= (1 - 9)(d{N) - 0d(N')) 

= (1 -9)dS e (N). 

We then come back to H and use this inequality to deduce that 


ce 


Ht 


ce 


Ht 


9 td{N) t - 9g' t d(N') t - —-— d(5gN) t = g+d(N) t - g^d(N) t - 9g' t d{N') t - —d(6 0 N) t 


ce 


\p\t 


< gtd6 e (N) t + 9(g+ - g' t )d{N') t - —d(S 0 N) t 


< gt dS e( N )t - 

< 0, 


7 


2(1 -9) 


d(S s N) t 


due to g + < g[ and g. < Hence dH t <C 0. To eliminate G, we set D t := exp (| J|t). By 
Ito’s formula, 

d{D t P t ) = D t [[P t J t - G t )dA t - dH t + cP t e ( S e Z t dM t + dS s N t )^j. 

But previous results show that (P t J t — Gt)dA t — dH t 0. Hence DP is a local submartin¬ 
gale. Thanks to the exponential moments integrability on \a\r and (Y')* (and hence | J|t), 
we use a localization procedure and easily deduce that 


Pt< E 


exp 


J s dA s ) Pt 




(23) 


We come back to the definition of Pt and observe that 


fc£<(i-0)|£|+«W£ 

<(i -0M- 


Hence (23) gives 


/ ~ e P\\M\ + \p\t \ 


< E 

< E 


exp ^ J J s dA s J exp {ce^ T 5 e () 
exp ( J J s dA s j exp ( 7 e 2/3||A|1 |£|) 


Ft 

F t 
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Hence 


1 _ Q 

SgY t < -InE 

7 


exp ^ 7 «e 2 ^ A II ^|£| 


(a s + 2/3\Y’\)dA ! 


Tt 


Therefore we obtain P-a.s. Y t < Y(, by sending 0 to 1. By the continuity of Y and Y', we 
also have P-a.s. Y. < Yf. 


□ 


As a byproduct, we can prove the existence of a unique solution given (A.3). 

Corollary 12 (Uniqueness) //(/,g,£) satisfies (A.3), P-a.s. g. > 0 and |£|, \q\t have 
exponential moments of all orders, then there exists a unique solution (Y, Z , N) to (/, g , £) 
such that Y* has exponential moments of all order and (Z ■ M T N) £ MP for all p > 1. 

Proof. The existence of a unique solution in the above sense is immediate from Theorem 
9 (existence), Theorem 11 (comparison theorem) and Corollary 10 (estimate). 


□ 


Remark. There are spaces to sharpen the uniqueness. The convexity in z motivates one 
to replace (A.3)(iv) by 

-QLt ~ 0\y\ - K\X t z\ < f(t,y,z) < a t + fi\y\ + ^\\t.z\ 2 . 

Secondly, in view of Delbaen et al [4] , we may prove uniqueness given weaker integrability, 
by characterizing the solution as the value process of a stochastic control problem. 


It turns out that a stability result also holds given convexity condition. The proof is a 
modification of Theorem 11 (comparison theorem). We set N° := N + U {0}. 

Proposition 13 (Stability) Let (/", g n , £™) ne N° with g n > 0 P-a.s. satisfy (A.3) associ¬ 
ated with (a", /3, 7 , ip), and (Y n , Z n , N n ) be their unique solutions in the sense of Corollary 
12, respectively. If - £° —► 0, f g T \/" - f°\(s, Y s °, Z° s )dA s —> 0 in probability, P-a.s. 
g n — 9° —t 0 as n goes to +oo and for each p > 0, 


sup E 

n6N» 


exp (p(|£”| + |a"|r)) 


sup | g n \ < 2 P-a.s. 
™eN° 1 


< +oo, 


(24) 


Then for each p > 1, 


limE 

n 


exp (p\Y n 



lim E 

n 



z° s ) T d(M) s (z: 


Z° s ) + d(N n 



Proof. By Corollary 10 (estimate), for any p > 1, 


sup E exp (p(Y")*) 
new 0 L 


Z:) T d{M) s Z^ + d(N n ). 


< Too. 


(25) 
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Hence the sequence of random variables 

exp(p|y"-y 0 r) +(J (( Z:-Z° s ) T d(M) s (Z:-Z° s ) + d(N n -N 0 ) s )y 

is uniformly integrable. Due to Vitali convergence, it is hence sufficient to prove that 
\Y n - y|* + J ((Z s ra - Z° s ) T d(M)(Z " - Z° s ) + d(N n - N) s ) —> 0 
in probability as n goes to +oo. 

(i). We prove u.c.p convergence of Y n — Y°. To this end we use 0-technique in the spirit 
of Theorem 11 (comparison theorem). For any 6 G (0,1), define 


8f? 

Sg n 

SgY n 


= f 0 (t,Y t °,Z?)-nt,Y t °,Zl 

= 9°-9 n , 

= y° - 9Y n , 




and 5gZ n ,5gN n , 5g(N) n , etc. analogously. Further, set 

f n {t,Y?,Z?)-f n (t,Y?,Z?) 


Pt ■— I { yO_yn^ 0 } : 
oP\\A\\ 


Y t ° - Y t " 


c := 




1-6 7 
P" := exp (ce Mt 6 e Y t n ), 

J- : = 7 eW« + 2^|yO|), 

D? := exp ( f J”^). 


Obviously p is bounded by /3 due to (A.3)(i). The 0-difference implies that 

f°(t,Y t °,z?)-er(t,Y t n ,z?) 

= Sf T + (0/"(t, y t °, zd - era, y t ", zd) + (/"(*,y t °, z°) - era,Yt . ^D)■ (26) 

By (A.3)(i)(ii), 

0 /"(t, y t °, ZD - ef n {t,Y t n , ZD = e Pt (y t ° - y t ") 


= Pt (0y t o -y t ° + y t o -0y”) 
<(i-0)^|y t o | + p t W, 

f n a, Y t °, Z t °) - 0/"(i, y t °, ZD < (1 - 0K + (1 - (9)/3|3f | + 




n 12 


2(1-0)'-^ 

Hence (26) gives 

/“(*, ^ t 0 , ^ t 0 ) - era, Y t n , zd - p t w < sr + a - «+ 2 ^ 1 ^°!) + 


i^zd 2 . 

(27) 
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To analyze the quadratic term concerning N° and N n , we deduce by the same arguments 
as in Theorem 11 that 


pp\p\t rP \p\t 

g?d(N°)t - Og?d(N n )t - —d(5 e N) t = 5g?d(N\ + g?d5 e (N)? - —d(6 e N n ) t 

« 9t (df>9(N)t - ^— e d(5 e N n )t) + 6g?d(N°) t 

< 6g?d(N°) t . (28) 

Given (27) and (28), we use an exponential transform which is analogous to that in Theorem 
11. This gives 


P t n < D?P? < E 




ye 


2/3||A|| r T 


1-9 


D:P?(\6f:\dA s + \6g:\d(N °) s ) 


Pt 


Using log a; < x and Y° — Y n < (1 — 9)\Y n \ + 8gY n , we deduce that 

D?P?(\8f?\dA a + \6g?\d{N°) t ) 


Y° - Y t n < (1 - 6>) \Y t n I + i—-E nn 


2/3||A|| r T 

D n pn . 

u t^t + l Q 


Tt 


Set 


A"(0) := exp (^^((F 0 )* + (F")*)) > P t ", 

/^ e 2/3||A|| 

S n (0) := exp (2__(| C o _ v ^ 


> Pp. 


We then have 


1 - 6 > 


F t °-F” < (l-0)|F t "| + ^E[i??H"(6l)+2^ :: ^D?A"(fl) J {\5f:\dA s +\8g:\d(N 0 ) s ) 
Now we use (A.3)(ii)(iii) to /" and proceed analogously to Theorem 11. This gives 


20\\A\\ 


1-9, 


7 


Y t n -Y t ° < (1-0)|F°| + —E D^(9)+^——D^A n (9) I (j \SadA s +\8g^\d(N°) a ) 


2/3||A|| 


l-6< 


Pt 


Pt 


Though looking symmetric, the two inequalities come from slightly different treatments for 
the 0-difference. The two estimates give 


1 - 9 


I Y t n - F t °| < (1 - 0)(|F t °| + 117*1) +-E £>?S"(0) 


Pt 


x; 


xf 


s 2 /3||^!Ie 


£>?A"(0) / (|<5/;|cL4 s + |<5 5 p|d(X 0 > s ) 


Pt 


A7 


We then prove w.c.p convergence of Y n — F u . For any e > 0, 


Y n _ y^oi* ^ 


e) < p((X x )* > |) + P((X 2 )* > |) + p((X 3 )* > |). (29) 
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We aim at showing that each term on the right-hand side of (29) converges to 0 if we send 
n to +oo first and then 9 to 1. To this end, we give some useful estimates. By Chebyshev’s 
inequality, 


p((x 1 )* > |) < ^■E[(y 0 )* + (w 1 )*], 

where E[(V 0 )* + (Y ra )*] is uniformly bounded. Secondly, Doob’s inequality yields 

P((X 2 )* > |) < . (30) 

Moreover, by Vitali convergence, the right-hand side of (30) satisfies 
limsupE< supE[(D”) 2 ]^ ■ limsupE[(S") 2 ]* 

n 

exp ^ 27 e 2 ^ j4 ^|£°|^ 

< +00. 


< su P E[(D") 2 ] 2 -E 


Hence, the first term and the second term on the right-hand side of (29) converge to 0 as n 
goes to +00 and 9 goes to 1. Finally, we claim that the third term on the right-hand side 
of (29) also converges. Indeed, Doob’s inequality and Holder’s inequality give 


D ((* 3 )* > f) 


< 


3e 2 /3||A|| 


< 


3 e 2 / 3 ||A|| 


E 


E 


D$A n (0) / (\8f:\dA s + \Sg:\d(N°) s ) 


(D^A n (9))‘ 


'E 


{\6f:\dA s + M\d(N°) ! 


2-i 


• (3D 


Note that 

f 

Jo 


(I Sf:\dA s + \6g:\d(N °) s ) < \cx\t + \a n \ T + 2||A||(F°)* + 7 (Z° • M + N°) T . 


Hence the left-hand side of this inequality has finite moments of all orders by Corollary 
10. Therefore, the left-hand side of (31) converges to 0 as n goes to +oo due to Vitali 
convergence. 

Finally, collecting these convergence results for each term in (29) gives the convergence 
of Y n -Y°. 

(ii). It remains to prove convergence of the martingale parts. By Ito’s formula, 


E 


L J o 


Z n s - Z° s ) T d(M) s {ZY - Z°) + d(N n - N°) s 

r T 


<E[|f*-£°| ] +2E \Y n -Y°\* / \F n {s,Y s n ,Z™) - F°(s,Y s °,Z° s )\dA a 


2E 


\Y n -Y°\* / ( g y(N n ) 3 -g 0 s d(N°) s ) 


As before, we conclude by Vitali convergence. 


□ 
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5 Change of Measure 


In the final section, we show that given exponential moments integrability, the martingale 
part Z-M+N, though not BMO, defines an equivalent change of measure, i.e., its stochastic 
exponential is a strictly positive martingale. We don’t require convexity which ensures 
uniqueness. But to derive the estimate for / Q T f(s, Y s , Z s )dA s , we use (A.2") where / is of 
linear growth in y. We keep assuming that P-a.s. \g. \ < The following result comes from 
Mocha and Westray [13]. 


Theorem 14 (Change of Measure) //(/, g,£) satisfies (A. 2") and |o|t have expo¬ 
nential moments of all orders, then for any solution ( Y , Z , N) such that Y has exponential 
moments of all orders and any \q\> £(q(Z ■ M + NJJ is a continuous martingale. 


Proof. We start by recalling Lemma 1.6. and Lemma 1.7., Kazamaki [9]: if M is a 
martingale such that 


sup E 
reT 


exp 


(vMr+(l 


-V )(M); 


< +oo, 


(32) 


for y ^ 1, then £[yM^ is a martingale. Moreover, if (32) holds for some y* > 1 then it 
holds for any y £ ( 1,77*). 

By Lemma 10 (estimate), Z-M + N is a continuous martingale. First of all, we apply 
the above criterion to M := q{Z ■ M + N ) for some fixed |<)| > Define A t (rj) such that 


In A t (??) := qy((Z ■ M) t + N t ) + q 2 (J^ - r^(Z ■ M + N) t . 


From the BSDE (2) and (A.2"), we obtain, for any r £ T, 

InA t ( 7 /) = qrj(Y t — Y 0 + J (f(s,Y s ,Z s )dA s + g s d(N) s )')+q 2 (J^-y)(Z-M + N) t 

< (2 + fi\\A\\)\q\ V Y* + |g|7?|a| T + \q\q{l + — ~ v)) {Z • M + N) T . (33) 


Note that 

7 , ]£l (1 
V 

Hence for any y > qo, (33) gives 


A T {y) < exp (\q\y(2 + /3)F* + \q\y\a\ T )- 


By exponential moments integrability, we have 

supE[A r (7y)] < + 00 . 
r£T 


It then follows from the first statement of the criterion that £ [qy{Z-M+N)') is a martingale 
for all y £ [go, oo)\{l}- The second statement ensures that it is a martingale for any y > 1. 
For any |g| > we set |g| £ (^, |g|), y := | > 1 , and apply the result above to conclude 
that £{q(yZ ■ M + N )) is a martingale. 


□ 
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